Clonal dynamics and copy number variants by single-cell analysis in leukemic evolution of myeloproliferative neoplasms.

Transformation from chronic (CP) to blast phase (BP) in myeloproliferative neoplasm (MPN) remains poorly characterized, and no specific mutation pattern has been highlighted. BP-MPN represents an unmet need, due to its refractoriness to treatment and dismal outcome. Taking advantage of the granularity provided by single-cell sequencing (SCS), we analyzed paired samples of CP and BP in 10 patients to map clonal trajectories and interrogate target copy number variants (CNVs). Already at diagnosis, MPN present as oligoclonal diseases with varying ratio of mutated and wild-type cells, including cases where normal hematopoiesis was entirely surmised by mutated clones. BP originated from increasing clonal complexity, either on top or independent of a driver mutation, through acquisition of novel mutations as well as accumulation of clones harboring multiple mutations, that were detected at CP by SCS but were missed by bulk sequencing. There were progressive copy-number imbalances from CP to BP, that configured distinct clonal profiles and identified recurrences in genes including NF1, TET2, and BCOR, suggesting an additional level of complexity and contribution to leukemic transformation. EZH2 emerged as the gene most frequently affected by single nucleotide and CNVs, that might result in EZH2/PRC2-mediated transcriptional deregulation, as supported by combined scATAC-seq and snRNA-seq analysis of the leukemic clone in a representative case. Overall, findings provided insights into the pathogenesis of MPN-BP, identified CNVs as a hitherto poorly characterized mechanism and point to EZH2 dysregulation as target. Serial assessment of clonal dynamics might potentially allow early detection of impending disease transformation, with therapeutic implications.

Inhibition of ERK1/2 signaling prevents bone marrow fibrosis by reducing osteopontin plasma levels in a myelofibrosis mouse model.

Clonal myeloproliferation and development of bone marrow (BM) fibrosis are the major pathogenetic events in myelofibrosis (MF). The identification of novel antifibrotic strategies is of utmost importance since the effectiveness of current therapies in reverting BM fibrosis is debated. We previously demonstrated that osteopontin (OPN) has a profibrotic role in MF by promoting mesenchymal stromal cells proliferation and collagen production. Moreover, increased plasma OPN correlated with higher BM fibrosis grade and inferior overall survival in MF patients. To understand whether OPN is a druggable target in MF, we assessed putative inhibitors of OPN expression in vitro and identified ERK1/2 as a major regulator of OPN production. Increased OPN plasma levels were associated with BM fibrosis development in the Romiplostim-induced MF mouse model. Moreover, ERK1/2 inhibition led to a remarkable reduction of OPN production and BM fibrosis in Romiplostim-treated mice. Strikingly, the antifibrotic effect of ERK1/2 inhibition can be mainly ascribed to the reduced OPN production since it could be recapitulated through the administration of anti-OPN neutralizing antibody. Our results demonstrate that OPN is a novel druggable target in MF and pave the way to antifibrotic therapies based on the inhibition of ERK1/2-driven OPN production or the neutralization of OPN activity.

Novel Molecular Insights into Leukemic Evolution of Myeloproliferative Neoplasms: A Single Cell Perspective.

Myeloproliferative neoplasms (MPNs) are clonal disorders originated by the serial acquisition of somatic mutations in hematopoietic stem/progenitor cells. The major clinical entities are represented by polycythemia vera (PV), essential thrombocythemia (ET), and primary myelofibrosis (PMF), that are caused by driver mutations affecting JAK2, MPL or CALR. Disease progression is related to molecular and clonal evolution. PV and ET can progress to secondary myelofibrosis (sMF) but can also evolve to secondary acute myeloid leukemia (sAML). PMF is associated with the highest frequency of leukemic transformation, which represents the main cause of death. sAML is associated with a dismal prognosis and clinical features that differ from those of de novo AML. The molecular landscape distinguishes sAML from de novo AML, since the most frequent hits involve TP53, epigenetic regulators, spliceosome modulators or signal transduction genes. Single cell genomic studies provide novel and accurate information about clonal architecture and mutation acquisition order, allowing the reconstruction of clonal dynamics and molecular events that accompany leukemic transformation. In this review, we examine our current understanding of the genomic heterogeneity in MPNs and how it affects disease progression and leukemic transformation. We focus on molecular events elicited by somatic mutations acquisition and discuss the emerging findings coming from single cell studies.

The Response to Oxidative Damage Correlates with Driver Mutations and Clinical Outcome in Patients with Myelofibrosis.

Myelofibrosis (MF) is the Philadelphia-negative myeloproliferative neoplasm characterized by the worst prognosis and no response to conventional therapy. Driver mutations in JAK2 and CALR impact on JAK-STAT pathway activation but also on the production of reactive oxygen species (ROS). ROS play a pivotal role in inflammation-induced oxidative damage to cellular components including DNA, therefore leading to greater genomic instability and promoting cell transformation. In order to unveil the role of driver mutations in oxidative stress, we assessed ROS levels in CD34+ hematopoietic stem/progenitor cells of MF patients. Our results demonstrated that ROS production in CD34+ cells from CALR-mutated MF patients is far greater compared with patients harboring JAK2 mutation, and this leads to increased oxidative DNA damage. Moreover, CALR-mutant cells show less superoxide dismutase (SOD) antioxidant activity than JAK2-mutated ones. Here, we show that high plasma levels of total antioxidant capacity (TAC) correlate with detrimental clinical features, such as high levels of lactate dehydrogenase (LDH) and circulating CD34+ cells. Moreover, in JAK2-mutated patients, high plasma level of TAC is also associated with a poor overall survival (OS), and multivariate analysis demonstrated that high TAC classification is an independent prognostic factor allowing the identification of patients with inferior OS in both DIPSS lowest and highest categories. Altogether, our data suggest that a different capability to respond to oxidative stress can be one of the mechanisms underlying disease progression of myelofibrosis.

Increased Plasma Levels of lncRNAs LINC01268, GAS5 and MALAT1 Correlate with Negative Prognostic Factors in Myelofibrosis.

Long non-coding RNAs (lncRNAs) have been recently described as key mediators in the development of hematological malignancies. In the last years, circulating lncRNAs have been proposed as a new class of non-invasive biomarkers for cancer diagnosis and prognosis and to predict treatment response. The present study is aimed to investigate the potential of circulating lncRNAs as non-invasive prognostic biomarkers in myelofibrosis (MF), the most severe among Philadelphia-negative myeloproliferative neoplasms. We detected increased levels of seven circulating lncRNAs in plasma samples of MF patients (n = 143), compared to healthy controls (n = 65). Among these, high levels of LINC01268, MALAT1 or GAS5 correlate with detrimental clinical variables, such as high count of leukocytes and CD34+ cells, severe grade of bone marrow fibrosis and presence of splenomegaly. Strikingly, high plasma levels of LINC01268 (p = 0.0018), GAS5 (p = 0.0008) or MALAT1 (p = 0.0348) are also associated with a poor overall-survival while high levels of LINC01268 correlate with a shorter leukemia-free-survival. Finally, multivariate analysis demonstrated that the plasma level of LINC01268 is an independent prognostic variable, suggesting that, if confirmed in future in an independent patients' cohort, it could be used for further studies to design an updated classification model for MF patients.

Gene expression profile correlates with molecular and clinical features in patients with myelofibrosis.

Myelofibrosis (MF) belongs to the family of classic Philadelphia-negative myeloproliferative neoplasms (MPNs). It can be primary myelofibrosis (PMF) or secondary myelofibrosis (SMF) evolving from polycythemia vera (PV) or essential thrombocythemia (ET). Despite the differences, PMF and SMF patients are currently managed in the same way, and prediction of survival is based on the same clinical and genetic features. In the last few years, interest has grown concerning the ability of gene expression profiles (GEPs) to provide valuable prognostic information. Here, we studied the GEPs of granulocytes from 114 patients with MF, using a microarray platform to identify correlations with patient characteristics and outcomes. Cox regression analysis led to the identification of 201 survival-related transcripts characterizing patients who are at high risk for death. High-risk patients identified by this gene signature displayed an inferior overall survival and leukemia-free survival, together with clinical and molecular detrimental features included in contemporary prognostic models, such as the presence of high molecular risk mutations. The high-risk group was enriched in post-PV and post-ET MF and JAK2V617F homozygous patients, whereas pre-PMF was more frequent in the low-risk group. These results demonstrate that GEPs in MF patients correlate with their molecular and clinical features, particularly their survival, and represent the proof of concept that GEPs might provide complementary prognostic information to be applied in clinical decision making.

Mutated clones driving leukemic transformation are already detectable at the single-cell level in CD34-positive cells in the chronic phase of primary myelofibrosis.

Disease progression of myeloproliferative neoplasms is the result of increased genomic complexity. Since the ability to predict disease evolution is crucial for clinical decisions, we studied single-cell genomics and transcriptomics of CD34-positive cells from a primary myelofibrosis (PMF) patient who progressed to acute myeloid leukemia (AML) while receiving Ruxolitinib. Single-cell genomics allowed the reconstruction of clonal hierarchy and demonstrated that TET2 was the first mutated gene while FLT3 was the last one. Disease evolution was accompanied by increased clonal heterogeneity and mutational rate, but clones carrying TP53 and FLT3 mutations were already present in the chronic phase. Single-cell transcriptomics unraveled repression of interferon signaling suggesting an immunosuppressive effect exerted by Ruxolitinib. Moreover, AML transformation was associated with a differentiative block and immune escape. These results suggest that single-cell analysis can unmask tumor heterogeneity and provide meaningful insights about PMF progression that might guide personalized therapy.

Characterization of New ATM Deletion Associated with Hereditary Breast Cancer.

Next-generation sequencing (NGS)-based cancer risk screening with multigene panels has become the most successful method for programming cancer prevention strategies. ATM germ-line heterozygosity has been described to increase tumor susceptibility. In particular, families carrying heterozygous germ-line variants of ATM gene have a 5- to 9-fold risk of developing breast cancer. Recent studies identified ATM as the second most mutated gene after CHEK2 in BRCA-negative patients. Nowadays, more than 170 missense variants and several truncating mutations have been identified in ATM gene. Here, we present the molecular characterization of a new ATM deletion, identified thanks to the CNV algorithm implemented in the NGS analysis pipeline. An automated workflow implementing the SOPHiA Genetics' Hereditary Cancer Solution (HCS) protocol was used to generate NGS libraries that were sequenced on Illumina MiSeq Platform. NGS data analysis allowed us to identify a new inactivating deletion of exons 19-27 of ATM gene. The deletion was characterized both at the DNA and RNA level.

Genomic Analysis of Hematopoietic Stem Cell at the Single-Cell Level: Optimization of Cell Fixation and Whole Genome Amplification (WGA) Protocol.

Single-cell genomics has become the method of choice for the study of heterogeneous cell populations and represents an elective application in defining the architecture and clonal evolution in hematological neoplasms. Reconstructing the clonal evolution of a neoplastic population therefore represents the main way to understand more deeply the pathogenesis of the neoplasm, but it is also a potential tool to understand the evolution of the tumor population with respect to its response to therapy. Pre-analytical phase for single-cell genomics analysis is crucial to obtain a cell population suitable for single-cell sorting, and whole genome amplification is required to obtain the necessary amount of DNA from a single cell in order to proceed with sequencing. Here, we evaluated the impact of different methods of cellular immunostaining, fixation and whole genome amplification on the efficiency and yield of single-cell sequencing.

Promoter Methylation Leads to Decreased ZFP36 Expression and Deregulated NLRP3 Inflammasome Activation in Psoriatic Fibroblasts.

The mRNA-destabilizing protein tristetraprolin (TTP), encoded by the ZFP36 gene, is known to be able to end inflammatory responses by directly targeting and destabilizing mRNAs encoding pro-inflammatory cytokines. We analyzed its role in psoriasis, a disease characterized by chronic inflammation. We observed that TTP is downregulated in fibroblasts deriving from psoriasis patients compared to those deriving from healthy individuals and that psoriatic fibroblasts exhibit abnormal inflammasome activity compared to their physiological counterpart. This phenomenon depends on TTP downregulation. In fact, following restoration, TTP is capable of directly targeting for degradation NLRP3 mRNA, thereby drastically decreasing inflammasome activation. Moreover, we provide evidence that ZFP36 undergoes methylation in psoriasis, by virtue of the presence of long stretches of CpG dinucleotides both in the promoter and the coding region. Besides confirming that a perturbation of TTP expression might underlie the pathogenesis of psoriasis, we suggest that deregulated inflammasome activity might play a role in the disease alongside deregulated cytokine expression.

Calreticulin Ins5 and Del52 mutations impair unfolded protein and oxidative stress responses in K562 cells expressing CALR mutants.

Somatic mutations of calreticulin (CALR) have been described in approximately 60-80% of JAK2 and MPL unmutated Essential Thrombocythemia and Primary Myelofibrosis patients. CALR is an endoplasmic reticulum (ER) chaperone responsible for proper protein folding and calcium retention. Recent data demonstrated that the TPO receptor (MPL) is essential for the development of CALR mutant-driven Myeloproliferative Neoplasms (MPNs). However, the precise mechanism of action of CALR mutants haven't been fully unraveled. In this study, we showed that CALR mutants impair the ability to respond to the ER stress and reduce the activation of the pro-apoptotic pathway of the unfolded protein response (UPR). Moreover, our data demonstrated that CALR mutations induce increased sensitivity to oxidative stress, leading to increase oxidative DNA damage. We finally demonstrated that the downmodulation of OXR1 in CALR-mutated cells could be one of the molecular mechanisms responsible for the increased sensitivity to oxidative stress mediated by mutant CALR. Altogether, our data identify novel mechanisms collaborating with MPL activation in CALR-mediated cellular transformation. CALR mutants negatively impact on the capability of cells to respond to oxidative stress leading to genomic instability and on the ability to react to ER stress, causing resistance to UPR-induced apoptosis.

Deregulated expression of miR-29a-3p, miR-494-3p and miR-660-5p affects sensitivity to tyrosine kinase inhibitors in CML leukemic stem cells.

The development of Imatinib mesylate (IM), which targets the oncogenic BCR-ABL fusion protein, has greatly improved the outcome of Chronic Myeloid Leukemia (CML) patients. However, BCR-ABL-positive progenitors can be detected in CML patients in complete cytogenetic response. Several evidence suggests that CML stem cells are intrinsically resistant to Tyrosine Kinase Inhibitors (TKI), and therefore they represent the most likely candidate responsible for disease relapse.In this work, we investigated the microRNA (miRNA) expression profile of different subpopulations of CML Leukemic Stem Cells (LSCs): Lin-CD34+CD38- and Lin-CD34-CD38- cells. These cell fractions have been previously shown to be endowed with TKI intrinsic resistance. Our analysis identified 33 common deregulated miRNAs in CML LSCs. Among those, 8 miRNAs were deregulated in CML independently from BCR-ABL kinase activity and therefore are likely to be involved in the BCR-ABL-independent resistance to TKI that characterizes CML LSCs. In particular, the up-regulation of miR-29a-3p and miR-660-5p observed in CML LSCs, led to the down-regulation of their respective targets TET2 and EPAS1 and conferred TKI-resistance to CML LSCs in vitro. On the other hand, miR-494-3p down-regulation in CML LSCs, leading to c-MYC up-regulation, was able to decrease TKI-induced apoptosis. These results demonstrate that aberrant miRNA expression in CML LSCs could contribute to the intrinsic TKI-resistance observed in these cell populations, and support the development of novel therapies aimed at targeting aberrantly regulated miRNAs or their targets in order to effectively eradicate CML LSCs.